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We propose a scenario for the prebiotic co-evolution of RNA and of fast folding proteins with large 
entropy gaps as observed today. We show from very general principles that the folding and unfolding 
of the proteins synthesized by RNA can function as a heat pump. Rock surfaces can facilitate the 
folding of amino acid chains having polar and hydrophobic residues, with an accompanying heat loss 
to the surrounding rock. These chains then absorb heat from the soup as they unfold. This opens 
the way to the enhancement of RNA replication rates, by the enzymatic action of folded proteins 
present in greater numbers at reduced temperatures. This gives an evolutionary advantage to those 
RNA coding amino acid sequences with non-degenerate folded states which would provide the most 
efficient refrigeration. 
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I. INTRODUCTION 

In order to understand the process of selection that 
could have lead to the observed properties of biolog- 
ical proteins, it is natural to look for mechanisms 
whereby proteins, which cannot duplicate themselves, 
evolved together with self-replicating molecules, "pre- 
living" concentrations of RNA, which provide the code for 
protein synthesis. [EI3iI3jI1iIMI3jm] The enzymatic action 
of proteins on the self-replication of RNA constitutes a 
hypercyle, which may be considered the elementary unit 
of evolution [l|, Q] ■ 

Biological proteins fulfill their functions in unique 
folded states, which they are able to reach in a very 
short time after being synthesized, [a] Small single do- 
main proteins fold into their secondary structures within 
milliseconds 9] , or even faster JJl] . These "native" states 
correspond to minimum (free) energy configurations. |ll| 
A random sequence of amino acids, however, will typi- 
cally have a degenerate ground state. (1,2, 13] We would 
like to propose a possible pathway in which those amino- 
acid sequences with unique ground states were selected 
in the course of evolution, without begging the question 
by referring to highly evolved biological functions. We 
will assume that this selection must have occurred in the 
prebiotic soup, where local temperature differences could 
have a large effect on the efficiency of RNA replication. 

As a chain folds (unfolds) at constant ambient tem- 
perature, heat will be given off (absorbed) by an amount 
proportional to the difference in entropy between the 
unfolded and folded states. Since we are interested in 
relatively high temperatures, and our proteins do not 
fold into their tertiary structures, we will ignore their 
hydrational entropy changes. For an amino acid chain 
of a given length, this difference will be the largest if 
the folded state corresponds to a non-degenerate ground 



state, or several possible low-energy configurations that 
are well isolated from each other by very high free energy 
barriers. We will argue that amino acid chains which es- 
sentially behave as two-state systems [J, ,14] with large 
entropy gaps prove to be the most effective refrigerants, 
if they were to be employed in a refrigeration cycle. Se- 
lection of amino-acid chains with folded states in deep 
free-energy wells could then be succeeded by the the evo- 
lution of more specific functions, leading to the pruning 
of those low-lying states so that only one, serving a highly 
specialized enzymatic activity, would survive. For con- 
venience, we will henceforth use the term "proteins" to 
mean amino acid chains, regardless of their degree of evo- 
lution. 

It has been pointed out [2] that a rudimentary form of 
compartition is necessary for evolutionary processes to be 
possible, and porous rock |ajl3 i^ among the likely envi- 
ronments to have played host to prebiotic processes, in a 
"soup," consisting of both organic and inorganic materi- 
als, that are assumed to be present in the prebiotic earth. 
For simplicity, let us concentrate on a series of compart- 
ments that contain RNA, protein molecules, amino acids 
and water. The surrounding rock is bathed by water at 
some medium temperature. The temperature range we 
have in mind is such that the RNA and proteins are sta- 
ble, and will vary, say, between 300-360 K, (the denatura- 
tion temperatures of most proteins are in this range j^])- 

In the next section. Section II, we will describe how 
proteins could act as a refrigerant in an adsorption re- 
frigerator. In Section III, we discuss a toy Hamiltonian 
for the guided and unguided folding of chains, and de- 
scribe the refrigeration cycle in the entropy-temperature 
plane. A discussion and pointers for further research are 
provided in the last section. 



II. PROTEINS IN AN ADSORPTION 
REFRIGERATION CYCLE 
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The hydrophobic interactions which drive the folding 
of proteins into their native states (Tj, Ha . HtI ] may also 
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FIG. 1: Cartoon showing the crossection of a pore, with the 
different steps in an adsorption-refrigeration cycle involving 
denatured and folded proteins. 



make it favorable for non-polar residues on the chain to 
adsorb on nearby hydrophobic surfaces jll| lowering the 
free energy of the whole system by reducing the num- 
ber of water molecules in interaction with the non-polar 
residues. If we assume that a certain fraction of the non- 
polar residues have adsorbed on the surface, in a rela- 
tively stretched conformation(as shown in Fig.l), then it 
can slightly increase the entropy with the diffusion and 
aggregation of the non-polar residues on the surface, al- 
lowing the intervening sections of the chain greater free- 
dom. The residues in the "loops" arc then free to fold 
into "beads" or "droplets" that are the incipient building 
blocks of the secondary structure. 

If this partially folded state is on the "correct path- 
way" 19, 20, 21] to a low energy folded conformation 
of the chain, stabilized by the specific intra-chain inter- 
actions for the given sequence of amino-acids, the rock 
surface can be said to act as a guide for the folding pro- 
cess. 

We would like to propose that the relatively large heat 
of adsorption and denaturation which can be achieved 
by proteins may make them amenable to functioning as 
refrigerants in an adsorption refrigerator. 

Unlike work-driven re frig erators, the conventional ad- 
sorption refrigerator 22, 23, 24, 25, 26| relics on the avail- 
ability of a cheap heat source and has no moving parts. 
The refrigerant absorbes heat from the cold reservoir at 
Tc, as it flash evaporates, and then is led to an adsorber 
bed (rather than a compressor), where it adsorbes on 
a substrate at some medium temperature T^, where it 
gives off a heat of adsorption Qads to the environment. To 
regenerate the refrigerant, the "low quality" heat source 
is used to heat the adsorber bed to a high temperature, 
T/j. The refrigerant that is released in this way is led 



to a condenser (which may be again at Tm), liquified, 
led through a nozzle to lower the boiling temperature, 
and then piped once more to the evaporation chamber. 
A "batch operating" system uses two "beds" (as adsor- 
ber and desorber), the cycle switching between them in 
turn. ^ 

The "gas" and "liquid" phases correspond, in our case, 
to the unfolded and folded states of the protein chains. 
The physical environment we have in mind is a series of 
interconnected chambers in porous rock bathed on the 
outside with ambient water. Part of the system is peri- 
odically exposed to intense heat. This may, for example 
be due to insolation or to the heat relased from a vent in 
the ocean floor. [2a | 

The refrigeration cycle consists of the following steps. 

i) The unfolded protein adsorbs on to a pore surface, 
where it folds. 

We assume the folding to take place isothermally, at 
the temperature of the rock surface, which we take to be 
r,„. An amount of heat Qads will be given off in this 
process. 

a) The regenerative step involves the heating of the 
surface to dislodge the folded proteins from the wall. 
This, we propose, may be supplied by sunshine or some 
geothermal source. During the regenerative step, a quan- 
tity of heat Qdcs will be absorbed. 

The chain can now detach from the rock, and com- 
plete its folding around its hydrophobic core, with the 
polar residues predominantly on the "outside." 0, y3 
We have shown [l^, using a lattice model for the hy- 
drophobic interactions, that above the temperature inter- 
val where the hydrophobic interactions make it favorable 
for a hydrophobic chain to adsorb onto a hydrophobic 
boundary, there is a temperature interval where the chain 
prefers to detach from the wall and go into a folded state. 
Raising the temperature further results in the unfolding 
of the chains. 

in) These proteins (denatured or otherwise) shall now 
be convected away from the hot wall. More temperate 
regions of the porous network, cooled by the ambient 
water, will act as the "condenser" in this system, and as 
they cool to the ambient temperature Tm once more, the 
proteins will equilibrate to the native, folded states. A 
quantity of heat Qcon is given off in the process. 

iv) In order for the denaturation temperature to be 
lowered as the proteins enter the cool compartment at Tc, 
(analogously to passing them through a nozzle and lower- 
ing the pressure, in a conventional refrigerator), here use 
is made of the fact that the denaturation temperature 
depends very sensitively on the total ionization strength. 
The gradient of ionization strength can, in fact, be easily 
achieved; the RNA molecules, the amino acids and dena- 
tured proteins in the interior of the chamber will lower 
the unfolding temperature for the coming chains. 

v) As the proteins unfold in the cool chamber due to 
the lowered denaturation temperature, they will absorb 
Qcv This phase of the cycle corresponds to the "evapo- 
rator" of the adsorption refrigerator. 



Clearly, for our adsorption refrigerator to work, the 
adsorber "beds" must be swept clean not just of proteins, 
but other unreacted amino acids, or other hydrophobes. 
The periodic nature of the heating serves this purpose as 
well. 

This refrigeration cycle is similar to magnetic cooling, 
if one thinks of the the action of the rock surface, facili- 
tating the folding, as the magnetic field. Note that for the 
chain to fold, and to give off heat to the medium tempera- 
ture reservoir, work has been done on it by the combined 
action of the rock surface (entropy mediated hydrophobic 
interactions) and intra-chain interactions. The ambient 
water bathing the outer surface of the rock acts as the 
thermal contact. 



III. A HIERARCHICAL TWO-STATE MODEL 
FOR PROTEIN FOLDING 

To investigate the refrigeration cycle in the entropy- 
temperature plane |27| . we used an exactly solvable toy 
system ^^ to model the temperature dependence of the 
entropy of the backbone of the protein chain as it folds 
or unfolds with or without guidance. Clearly this is very 
schematic picture, but we believe it conveys the essential 
physics. 

We assume that there is only one folding pathway lead- 
ing to a single low-free energy folded state, to simplify the 
discussion. In this respect, then, the amino-acid chains 
under consideration are similar to modern day proteins. 

The existence of a unique pathway means that an or- 
dered sequence of binding events occur between different 
parts of the protein [2l| ; and if this particular sequence 
is not followed, the protein can not fold. 

Slightly modifying the hierarchical zipper-like model 
proposed by Bakk et al. |23|, let us assign a variable 
CTj = 1, ... g to the different choices that can be made at 
each node, with only say a* leading to a correct folding 
move. Then, the state variable at the ith node on the 
folding pathway of N nodes [20| may be written in terms 
of a Kroenecker delta as, 



The Hamiltonian is 

N 



where 



7^ = -Ae ^ *, - (1 - A)e*Ar 



*» = n ^fc 



(1) 



(2) 



(3) 



For A 7^ we see that intermediate partially folded 
states also lower the energy, as would be the case under 
the guidance of chaperons [3 j by an amount Ae while the 
energy gap of the native state is given by e. For A = 0, 
this Hamiltonian allows the protein to be in two distinct 



states only, native and unfolded. No unfolding can occur 
inside an already folded part of the protein. Notice that 
(PI differs from the Bakk et al. model in the last term, 
which in our case is not multiplied by N] this allows the 
folding transition to shift to higher temperatures under 
guidence, as is experimentally observed p33 |. 

The partition function for this model can be evaluated 
exactly. For a protein having a folding pathway of N 
nodes, with (3 = l/fcsT, ks the Boltzmann constant, 
and T the temperature. 



Z = iq-1 



)qA'-i (^j^) + e^[^(^-i)+il^ . (4) 



where we have defined ly = (3Xe — \nq. (Note that 
L'Hopital's rule must be used in the first term in case 
the denominator vanishes). From the Hclmholtz free 
energy F = — kT In Z one may compute the entropy 
S = -dF/dT)v and heat capacity C = TdS/dT)v. (In 
this model the free energy is independent of the volume, 
therefore the partial derivatives above can be treated as 
ordinary derivatives.) The entropy is plotted in Fig. 2 
for different values of A, A = corresponding to the 
unguided and A 7^ to the guided cases. The system 
exhibits a sharp transition (becoming first-order in the 
thermodynamic limit, i.e. limA'^ — > 00) for the unguided 
case, as expected from this two- state model. The effect 
of guiding is to lower the entropy, smooth the transition 
and to shift it to higher temperatures, and is similar to 
the effect of turning on a field, in the case of magnetic 
phase transitions. 

The refrigeration cycle for the refrigerant is shown 
in Fig. 2 in the entropy-temperature plane. The sharp 
unfolding transition (1 -^ 1') analogous to the flash- 
evaporation in a conventional adsorption refrigerator, 
takes place in the cooling chamber, at Tc- The high pH 
conditions we take to correspond to the A = curve. 
The denatured proteins are convected towards the ad- 
sorbing walls (the walls could be coated by some lipids 
to make them more hydrophobic 28i|), warming up in the 
process(l' — ^ 2), and adsorb (2^3). Upon adsorption A 
is set to unity, and the entropy of the chain drops. This 
happens at the (fixed) intermediate temperature Tm of 
the adsorbing wall, and the heat of adsorption is carried 
away by the water bathing the outer wall of the compart- 
ment. When the heat is turned on, as it periodically is, 
the wall will heat up to T^, along the curve 3^4. In 
the process, Qdes is absorbed by the proteins. At point 
4 in the cycle, they become free of the wall; and are now 
convected along 4 ^- 4' ^- 1, back to the starting point. 
This happens along a curve with a nonzero value of A, 
here chosen to be A = 0.13, because, as explained in 
item (ii) of the previous section, we assume these pro- 
teins remain near the hydrophobic wall as they give off 
excess heat to the rock boundary washed by cooling wa- 
ter on the outside, and are convected back towards the 
cool chamber. 

The various heats of evaporation, adsorption, desorp- 
tion and condensation may, in principle, be computed 




FIG. 2: The refrigeration cycle in the entropy-temperature 
plane. The curves have been drawn for A = 0, 0.13 and 1. 
For this plot q — 8 and A'^ = 10. The temperature axis is in 
units of e/fcs. 



from the entropy as a function of the temperature for 
the various values of A. It is clear from the geometry 
of the curves in Fig. 2 that the area to the left of the 
curve between any two points (ij) on the curve may be 
obtained by doing the integral over T instead, viz., 



for typical proteins, where L, the length of a protein 
chain, is taken to be 100, T^ = 300 K and fc = IQ-^/s. 

On the other hand typical heat conductivities for rock 
are around Kt = 2 - 6 J/m-K-s jM EJ. With a tem- 
perature difference between the ambient water and the 
inside of the cooling chamber being, say T^ — Tc ~ 30 
K, and a surface area ~ r^, the heat conducted into the 
cooling chamber from the surroundings, Jin — dQm/dt 
will be 



Jin = Kt 



T„ 



(8) 



where £ is the thickness of the wall seperating the cham- 
ber from the ambient water. With r^ ^ \ cm, and taking 
JsTt = 4 J/m-K-s, one finds Ji„ - (1.20/£)J/s, if (. is ex- 
pressed in cm. . The break even point where Jjn — Jov 
is around £ '^ 2.5 cm. 

The power needed for this system can be found from 
Jdcs — dQdcs/dt. The insolation at the earth surface is 
of the order of Wso\ — 200 J/m^-s. Requiring 



Js = VKsol^b > Jdcs = 



Jr.' 



COP 



(9) 



where rb is the radius of the adsorption bed receiving the 
solar power and inserting numbers, one finds rb should 
be of the order of a few tens of cm, for COP as low as 
0.1. 



y«j ^ ^j'^x 



T,S 



(i) 



SxiT)dT 



(5) 



The "coefficient of performance" of the refrigerator is 
given by |2^ 



COP = 



Qdc 



Wo 



(6) 



where Wg is a small amount of work normally needed 
to operate a pump. In our case, this is provided by the 
gravitational force driving the convection current. To 
maximize Qcv, one must have essentially two-state un- 
folding at Tc, with a maximal entropy gap. On the other 
hand, one must clearly have the folding-unfolding rates, 
Qcon ~ Qcv < Qdcs — Qads to rcspcct the first law of ther- 
modynamics. For this, one needs the adsorption to take 
place at r,„ values that are as low as possible, while Qcon 
should also be minimized. We have chosen appropriate 
values of A in Fig. 2 to make this plausable. 

In order to have a crude estimate of the actual pump- 
ing rates possible, one should note that the entropy of de- 
naturation of a typical protein is about AS' ~17J/mol-K 
per residue. |33,|3J| The folding-unfolding rates range be- 
tween k ^ 10^/s to 10^/s '10, 35] whereas typical protein 
concentrations, namely p, are of the order 10/iM lOj. If 
one takes a cooling chamber of linear dimensions Tc ~ 1 
cm, one finds that the heat pumped away, Jov — dQcv/dt 
is 



Jov = pkrlLTcAS - 5 X lO^^J/s 



(7) 



IV. DISCUSSION 

We have shown above that a protein soup within a 
porous rock could function as a self-regulatory refrigera- 
tion cycle and lower the temperature of the soup within 
the pore, for realistic ranges of the physical and chemical 
parameters. The efficiency of the cycle strongly depends 
on the size of the entropy gaps. 

It is well known that the presence of protein molecules 
acting as enzymes may effect the RNA replication rates 
(which we have not considered) by factors of up to 
10^ 0, [33. Catalytic activity is a function of the spa- 
tial structure [l| and therefore requires the proteins to 
be in a unique folded state, whose stability also depends 
on the temperature and is optimized only for a definite 
temperature interval [3,|33. Moreover, RNA replication 
rates depend non-monotonically on the temperature |38l | 
and drop off outside a definite temperatre range. We 
have shown that proteins with large entropy gaps are able 
to achive temperatures in a small compartment that are 
lowered relative to the ambient temperature. This vindi- 
cates our initial assumption that this criterion could con- 
cievably have played a role in their selection, in a environ- 
ments that are too warm for the optimal self-replication 
of RNA. 

There is experimental evidence that the folding tran- 
sition is like a two-state system for many single-domain 
proteins |9l ll4 . l3]| . Although the proteins fold via a two- 
state folding pathway, especially at higher temperatures. 



the presence of some intermediate states might be neces- 
sary for the folding process to find the native state. It has 
been experimentally observed [22| that the concentration 
of chaperons in E. coli rises as the temperature increases, 
indicating that at high temperatures E. coli needs help 
in order to fold its proteins. 

The efficiency of protein folding can be adversely af- 
fected if partially folded proteins aggregate in order to re- 
duce exposed hydrophobic residues. Molecular chaperons 
bind reversibly to these partially folded chains prevent- 
ing their aggregation and promoting their passage down 
the folding pathway. Therefore, one might speculate that 
these chaperons have taken the role of rock surfaces in 
the course of higher evolution. 

It is interesting to note that there exist so called heat- 
shock proteins which are synthesized in large numbers 
when the temperature of the environment is suddenly 
raised above (or dropped below) those temperature above 
and below which proteins in the cells would normally de- 
nature, and chaperone the correct folding of other pro- 
teins in the cell. [32| These proteins observed today under 



extremely hot conditions, in archeabacteria are relatively 
small and have very fast folding rates. 0,l3,ll3 These ex- 
amples seem to point also in the direction of a direct cor- 
relation between thermal properties of the environment 
and protein functions, which may have evolved from the 
rudimentary function proposed here. 

Other studies also have demonstrated that rock sur- 
faces may have played a selection role in prebiotic con- 
ditions, c.f., the selection of right handed amino acids 
binding to optically active surfaces of calcite crystals with 
same chirality |40j. 

Clearly we are making no claims that the present pro- 
teins are identical to the end products of one particular 
selection mechanism. Once the temperature was suffi- 
ciently lowered, other, more complex evolutionary pres- 
sures would come into play. However, our calculations 
show that the proposed refrigeration cycle could very well 
have played an important role in the co-evolution of fast 
folding proteins with large entropy gaps and the RNA 
molecules which code them. 
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